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Abstract: We introduce a family of matrices with non-commutative entries that generalize the classical real Wishart 
matrices. With the help of the Brauer product, we derive a non-asymptotic expression for the moments of traces of 
monomials in such matrices; the expression is quite similar to the formula derived in [51 Theorem 2.1] for independent 
complex Wishart matrices. We then analyze the fluctuations about the Marchenko-Pastur law. We show that after 
QQ ' centering by the mean, traces of real symmetric polynomials in g- Wishart matrices converge in distribution, and we 
' identify the asymptotic law as the normal law when g = 1, and as the semicircle law when g = 0. 
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1 Introduction 



The real Wishart distribution was introduced by Wishart [46]. As evidenced by the vast Hterature, the Wishart 
law is of primary importance to statistics, see e.g. [Tl[3S]. The compound Wishart distribution was introduced 
by Spcicher |41j . We will consider a slightly more general class of laws on random matrices which was introduced 
in Rcf. [20j . and also its non-commutative generalization. 

Definition 1.1. Let S be an N x N real positive- definite matrix (i.e. x*Sx > for all non-zero x £ M.^ ). Let 
^ , B be an M X M real matrix. We will say that a random matrix W is Wishart with shape parameter B and scale 
' parameter S, if 

W^A'X'BXA, (1.1) 

, where X is an M x N matrix of i.i.d. N{0, 1) random variables, and A is a symmetric root. A' A = A^ = T,. 
■ (Here and throughout the paper, A' denotes the transpose of A.) We write W G W(S,i?). 

o ; 

OO . When B is positive-definite, Ref. (TU] interprets W as a sample covariance under correlated sampling. 

Similar interpretation for the complex case appears in Ref. [TT]. Somewhat more generally, one could consider 
all matrices of the form (|l.ip with non-symmetric A; the moments of such matrices would then depend on 
E+ = A' A and = AA' as well as on B, B' . 

The usual real Wishart matrices correspond to the choice B ~ Lm, so the notation W{'E,M) in [U Section 
[ 7.2] means >V(S,/a,/) in our notation. Ref. [17] uses a slightly different parametrization a = 21], p = M/2, 
which should be kept in mind when comparing the formulas. Compound Wishart matrices are the elements of 
W{In,B), see [HI Section 4.4] or [HI page 169]. 

Our main result about real Wishart matrices. Theorem 12.91 gives a non-asymptotic expression for the 
moments of traces of monomials in several independent real Wishart matrices. This formula is similar to the 
expression derived in [9] Theorem 2.1] for independent complex Wishart matrices. The statement requires 
additional notation on Brauer products which will be introduced in Section [2l Our main application of this 
formula deals with the fluctuations about the Marchenko-Pastur law. The latter problem fits naturally into a 
more general non-commutative setting, see Theorem 11.71 

We now introduce matrices which arise from matrices with noncommutativc g-Gaussian entries by a formula 
similar to (jl.l|) . The noncommutativc (/-Gaussian random variables were introduced on the algebraic level by 
Frisch and Bourret [l^ ; additional physical applications are in [38l [44] . The Fock-space construction is due to 
Bozejko and Spcicher [7]; an excellent exposition of this topic appears in Ref. [6j. 

Fix a parameter — 1 < <? < 1. Let ,4 be a noncommutativc probability space, i.e., a unital complex *-algebra 
with a tracial state t : A ^ C; see, e. g., [U] Sect. 1.2] or [35] for more details and examples. 
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Definition 1.2 ( [Hi Theorem 1]). Xi, X2, ■ ■ ■ (z A are q-Gaussian, ifK-j = X*, and for any finite collection of 
Yi, . . . , Y„ e {Xi, X2, . . . }, 

The sum is over all pair partitions V of {1, . . . ,n = 2k} and to be interpreted as when n is odd; cr(V) := 
^2) ■ {Hi ji}, {^2, ^2} S V, ii < Z2 < ji < ^2} is the number of crossings of a pair partition V. 

Consider now the *-probabihty space Mnxn{A) of ah N x N matrices with elements from A, with matrix 
multiplication as the product operation, with state 

1 ^ 

r o trAT (X) ■■= j^Yl ^( 

i=l 

and with the natural definition of the star operation as conjugation of elements in A, and transposition. 

For > 1, a positive definite real N x N matrix S and an AI x M positive definite real matrix B, let X^.j 
be g-Gaussian noncommutative random variables such that for 1 < i, A: < M , 1 < j,m < N ^ 

T(X,jXfe,,„) = (1.2) 

Definition 1.3. We will say that W G A4nxn{A) is a q-Wishart random variable with shape parameter B and 
scale parameter S if 

W = X*X, (1.3) 

where X = [X^j] is the M x N matrix with q-Gaussian entries that satisfy p.2p . We will write W G Wg (S, B). 

(Matrices with noncommutative entries have been studied in [Ml [351 [Mj ■ According to [301 Definition 

I. 2.2] when E = /jv, -B = /m, X^j- form a g-circular system.) 

We will also need a substitute for independence. 

Definition 1.4. We will say that Wi, W2, . . . , G Wg (S, S) are q-orthogonal if the entries of the 
corresponding matrices (|1.3p ; {X^j^j} for Wi, {X^^^^} /or W2, are pairwise uncorrelated, T(X|^^X^fy) = 

Remark 1.5. When q = 1, all g-Gaussian random variables commute and their joint moments are the moments of 
the multivariate normal law of matching covariance. In particular, for every finite q-Gaussian family Xi, . . . , X^, 
there exists an m x m real matrix A and independent classical standard normal random variables Zi, . . . , Zm 
such that with Yj = J^ki^hk^k, 

t(X^ . . . Xt) = ]E(r^ . . . Y^r) for aU fci, . . . , fc^ > 0. 

Thus the joint moments of g-orthogonal Wi, W2, . . . , G Wg (S, S) coincide with the respective joint 
moments of independent real Wishart matrices. Under appropriate assumptions, such matrices are asymp- 
totically free and their limit laws are compound Marchenko Pastur laws. The conclusion of Theorem 1 1.71 in this 
case means that for real Wishart matrices the fluctuations around the Marchenko-Pastur law are asymptotically 
normal. (This fact is known for complex Wishart matrices, see |33j.) 

Remark 1.6. When q = 0, for every finite q-Gaussian family Xi, . . . , X,„, there exists an m x m real matrix A 
and free semicircular elements Zi, . . . , Z„j such that with Yj = '^i;[A]jk'Zk, 

t{Xj, . . . X,„ ) - r (Yj, . . . Y,„ ) for ah 1 < , . . . , j„ < m. 

Thus, for example, the traces of powers of q-orthogonal Wi, W2, . . . , G yVg(S,i?) are free in {A,t). In 
Proposition 12.191 we show that W G Wg (-^/at, B) as an element of {Mnxn{A),t o Hn) has a compound 
Marchenko-Pastur law already for finite N, so on average the fluctuations around the limit are zero. Theorem 

II. 71 shows that the un-averaged fluctuations around the Marchenko-Pastur law have a non-trivial semicircular 
limit. 

Asymptotic normality for traces of polynomials in several independent complex Wishart matrices is 
established in Ref. [33] . This case served as a departure point for a number of papers on second order freeness, 
see e.g. [12] and the references therein. For an extension to more general covariances see also [H Theorem 1]. 
Asymptotic normality for one real Wishart matrix is established in ^ Theorem 2], [531 Theorem 4.1]. Asymptotic 
normality of self-adjoint polynomials in noncommutative entries in a different setting appears in [241 Theorem 
1]. It turns out that in our noncommutative setting, asymptotic limit law always exists, but it depends on the 
polynomial in a nontrivial way. 



Theorem 1.7. Suppose Wi^\ W^^\ . . . , e ( jjIn^Im) are q-orthogonal and M ~ M (N) is such that 

M/N ^ A G [0, oo) and N ^ oo. If Q{) is a fixed non- commutative polynomial in s variables, and 



then Xn converges in moments as — > oo. Furthermore, if Q is a real symmetric polynomial so that is 
self-adjoint, then Xm converges in distribution; if in addition q = 1 then Xn has asymptotically normal law; if 
q ~ then X^ has asymptotically semicircle law. 

Theorem 11.71 raises an interesting question about how to identify or describe the hmiting law for aU 
q Cz [—1,1]. Our proof rehes on a general property of (/-Wishart matrices, see Theorem 14.11 and freeness or 
independence when g = or g = 1; this technique is not helpful in general. 

Computer-assisted calculations based on exact formulas from Theorem 12.181 indicate that the limit law for 
general q depends on the polynomial in a complicated way, see Example 14.81 Similar situation was noted for 
Wigner matrices in noncommutative entries in an unpublished manuscript [31j . 

For one matrix case, one may be able to find polynomials with limits that can still be described in simple 
terms. One such polynomial seems to be tr(W'-^^), see Example 14.61 another suggestive computer-assisted 
calculation is in Example 14.71 

We now list two corollaries. We begin with traces of polynomials of a single g-Wishart matrix. 

Corollary 1.8. If Q is a real polynomial, W G Wq [jjIn,Im) , M/N X E [0, oo), then 



converges in distribution as N ~> oo. Furthermore, if q = 1 then X^ has asymptotically normal law; if q = 
then Xn has asymptotically semicircle law. 

For 5 = 1, the traces commute, and their joint moments coincide with the moments of traces of the 
corresponding real Wishart matrices (HI]); taking QiW^, . . . , WJ = Y^^i n'=i ^t(r), by the Cramer- Wald 
device ([H Theorem 29.4]) we deduce the following multimatrix generalization of Theorem 2], [531 Theorem 
4.1] which shows that [331 Corollary 9.4] holds also for real Wishart matrices. 

Corollary 1.9. Fix s > 1 and t -.N {1, . . . , s}. Suppose Wi, . . . , G W Im) are independent. For 

a fixed k > I that does not depend on N , and for 1 < j < k consider real random variables 



converges in distribution to the multivariate normal law. 

The paper is organized as follows. In Section [5] we introduce notation, state our results on mixed moments 
of traces of monomials in independent Wishart matrices and on mixed moments of traces of monomials in 
g-orthogonal q- Wishart matrices, and we deduce some of the related expressions available in the literature. 
Section [3] contains proofs of the formulas for moments. Section |3] contains some auxiliary results, and the proof 
of Theorem 11.71 

2 Moments of real Wishart and g- Wishart matrices 

Explicit formulas for moments of polynomials in one real Wishart matrix appear in [T71 dHl HOI HE] • The formulas 
are more complicated than the formulas for the complex case [Ml [20l [33] , and often involve sophisticated tools, 
like wreath products. Jack polynomials, and properties of hyperoctahedral group. 

One of the main results of this paper is the non-asymptotic expression for moments of traces of monomials 
in several independent real Wishart matrices, and its noncommutative q- generalization. In Section [H we use the 
moment formula to prove Theorem 11.71 but since the formula is of independent interest we consider a more 
general case than what we need for the proof. Our approach follows closely [9l Theorem 2]: it relies on Wick 
formula, combinatorics from [33] and the Brauer product [B] of pair partitions. The connection of the Brauer 



X^ -.^ tr(Q(Wr\ . . . , WW)) - r (tr(Q(w(^), . . . , W^))) , 



XAr:=tr(Q(W))-T(tr(Q(W))) 




Then as N —^ oo, M/N ^ A G [0, oo), the sequence of k- dimensional random variables 




algebra to integration over the orthogonal group has already been noted in il-il who seem not to use the Brauer 
product explicitly. With the use of the Brauer product, our formula closely resembles the formula for the complex 
case. 

Following the referee's comment we remark that we use Brauer algebra as a convenient notational 
replacement for the permutation group, but since we multiply only by elements from our use of the Brauer 
structure is in fact the action of Sn on 2n-pair partitions, as in [17| . In particular, the constant of the Brauer 
algebra does not enter our formulas. 

After the first version of this paper was written, we learned from J. Mingo about his work [32] with Emily 
Redelmeier who has found a very elegant diagrammatic interpretation of the right hand side of equation ()2.9p . 

2.1 Pair partitions, permutations, Brauer product 

Consider matrix variables xi, X2, . . . , Xs which we interpret as variables of "colors" 1, . . . , s. We are interested 
in monomials p^,t(xi, X2, . . . , x^) of degree n which we parameterize by perfect matchings (pair partitions) 7 of 
{±1, ±2, . . . , ±n}, and by functions t : {1, . . . , 71} {1, . . . , s}, which we interpret as assigning "colors" to the 
integers. To introduce these polynomials in formula (|2.6p below we need additional notation and terminology. 

2.1.1 Pair partitions 

Let J-n denote the set of pair partitions (i.e., partitions into two-element sets) of the 2n-clcment set 
{±1, ±2, . . . , ±n}. Ref. [TS] uses the notation {1, 1, 2, 2, . . . , n, n} for the same object - the ordered set of 2n 
elements with the distinguished match. Other authors interpret J-"„ as the set of 1-factors, or perfect matchings, 
or 1-regular graphs on the vertices {±1, ±2, . . . , ±n}. 

Identifying the pairs of a pair partition 7 6 with the cycles of the permutation, we can embed 

J^n C <5{-|-i -|-„} = S2n- Here Sp denotes the group of permutations (i.e., bijections with composition) of a finite 
set F. As usual 5„ = 5{i In the embedding, J^n is the set of all involutions on {±1, . . . , ±n} with no fixed 

points. Accordingly, we write 7(1) = j if {i, j} £ 7. 

We will also find it convenient to represent 7 as a graph with vertices arranged in two rows 

1 2 ... n 

-1 -2 ... -n 

with the edges drawn between the vertices in each pair of partition 7. For example, a pair partition 5 which 
pairs j with — j for j = 1, . . . , n is identified with the permutation 5{j) — — j, or with the graph 



(2.1) 

In fact, S is the distinguished matching on {±1, . . . , ±n} and plays a special role in several definitions that 
follow. 

With each 7 G JF„ we associate the permutation 7r(7) G iS„ and the sequence £(7) of ±. These mappings 
TT : J-'n — > Sn and e : J-n { — , +}" are defined in terms of graphs as follows. 

Consider the 2-regular graph 5 U 7. To define the first cycle of permutation a :~ 71(7), start with vertex — 1. 
Follow upwards the first edge of S towards 1 , and then continue by following the consecutive edges of the graph 
(5 U 7. The first cycle of a consists of the vertices from the upper row {1, 2, . . . , n} that were visited by this path, 
taken in the order in which they were visited. The corresponding values of £(7) on this cycle are determined by 
the direction of traversing the edges of S: we assign + if the corresponding edge of S was traversed upwards, and 
— otherwise. (In particular, the sequence £(7) of signs always starts with +.) 

The next cycle of 71(7) consists of the entries in the first row, starting with the leftmost element j of the 
first row that was has not yet been visited. To determine the cycle, we again start at the corresponding vertex 
—j from the lower row, follow upwards the edge of S towards j, and continue in the same direction along the 
edges of 5 U 7 until we complete the second cycle. We then repeat these steps until we exhaust all the elements 
of the first row. (Similar construction appears in Refs. [TKl [T71 [^ .) 

For example, suppose 




7 = {{l,2},{3,-4},{-l,-2},{-3,4}}= ' ' 



Then 



tx 



Here the edges for the initial vertices of each cycle are oriented in the direction of the "first step". Thus 
7r(7) = (1,2)(3, 4) and £(7) = (+,-,+,+). 

Next, we introduce partitions that will play a role in Section [ 



(2.2) 



Example 2.1. For A h n, i.e. for X ~ (Ai > A2 > • • • > 0) with ^ Xj ~ n, consider 

{l,-2},{2,-3},...,{Ai,-l} , {Ai + l,-Ai-2},...,{Ai + A2,-Ai-l} 



Ai pairs 



A2 pairs 



Then 




Ai + l Ai+2 



A1+A2 A1+A2 + I 



and 



Ai+2 




A1+A2 



A1+A2 + I 

t 



Thus tt{<7\) ~ (1, 2, . . . , Ai)(Ai + 1, . . . , Ai + A2), . . . has cycle type X. In particular, 6 — ai^ and 7r((5) = id. 
Similarly, 7r(cr„) ~ (1, 2, . . . , n). 

Definition 2.2. By J-^ we denote the set of pair partitions a G -F„ in which the vertices of the upper row 
are always connected to vertices of the lower row. Such pair partitions are characterized by the fact that 
e=i+, +,...,+). 

Remark 2.3. Brauer [8] identifies the elements of J-'^ with permutations, as tt : J-^ ^ <S„ is a bijcction. For 
a G the definition of Tr{a) simplifies to 

7r(a)(i) -.^-aij), j-l,2,...,n. (2.3) 

The inverse mapping iS„ 9 a 1-^ cr G is then a{j) — ~a{j). 

2.1.2 Brauer product of pair partitions 

The Brauer product of pair partitions has simple description in terms of graphs. To define cr ® 7 we draw the 
two rows for 7, then we draw the graph for a so that its upper row covers the lower row of 7. We construct 
the new two-row graph for cr © 7 by removing the middle row, but retaining the edges between the vertices of 
the first and third row that were connected in the three-row graph. Brauer [8] weights the resulting graph by a 
coefficient that depends on the number of cycles lost in the middle row; we will apply the Brauer product only 
when there are no such cycles, so our weight will always be 1. 

Example 2.4. Suppose 



a = {{1,-2}, {2, -3}, {3, -4}, {4,-1}} = 




7 = {{1, 2}, {3, -4}, {-1, -2}, {-3, 4}} = 
The intermediate three-row graph is 

• — • • • 

X 



This gives 

• — • • • 

^©7= I ={{l,2},{-2,-3},{-l,3},{-4,4}}. 

• • — • • 

Our formulas need only the Brauer products with cr G J-^ , and in this case the analytical definition takes 
the following form. 

Definition 2.5. For a G J-^ and 7 e JF„, the Brauer product 0/ a @ "f is defined as an element of Tn such that 

'<y{-i{i)) */.? > 0, 7(j) < 0, 
^ '^^ ^ <y{-i{-<ym < 0, 7(-^(j)) < 0, ^ ' 

if 3 < 0, 7(-^(j)) > 0. 

(These are the same cases that are listed in Table [2j) 
Remark 2.6 (Brauer). If (Ti, (T2 G J-n then tri ® (T2 G J-^ and 

7r((Ti ® 0-2) = 7r(cri) o 7r((72) (2.5) 

becomes just the ordinary composition of permutations. 

We also remark that ^ ® 7 = 7 and, when the more general definition is used, 7 ® (5 = 7. 
It will be convenient to have concise notation for certain objects and properties of 7 G JF„ which depend 
only through the associated permutation 7r(7) G Sn- 

Definition 2.7. For a "coloring" t : {1, . . . ,n} ^ {1, . . . s}, by 

T„{t) = {7 e :f„ : tMj)\) = t{\j\),j = ±1, . . . , ±n} 

we denote the set of color-preserving pair partitions. Such partitions are characterized by the property that the 
permutation a :— 7r(7) is color-preserving, t o a = t, see JEt (4)]- 

By C{j) wc denote the set of cycles of permutation 7r(7). For 7 G J-n{t), all cycles have the same color; by 
Cj(7) we denote the set of cycles of irij) which are of color j G {1, . . . , s}. 

For a 7 G .?>i(t), we can color each edge of 7 by the color of the adjacent vertices. Since for a G T:^ each 
edge of cr ® 7 arises from combining exactly one edge of 7 with some of the edges of ct, we can color each edge 
of CT ® 7 by the color inherited from this unique edge of 7 G !Fn{t)- 

Definition 2.8. For j £ J-n{t) , a £ J-^ , byt{a,^) we denotethe coloring of {1,2, ... obtained as follows. For 
1 < i < n, we define t{a, 7)(i) as the color of the edge of a © j which is crossed on (5 U (cr ® 7) while connecting 
i to the next element of the corresponding cycle of it (a @ 7). 

2.1.3 Formula for moments 

For a "coloring" ..,??.} ^{l,...s}, a pair partition 7 G J-n, and matrix- variables Xi, X2, . . . , Xg, define 

monomial 

p^,t(xi,X2,...,x,s) = Jl ti' n^Hj) ■ (2-6) 

cec(7) / 

Here the notation tr ^Iljec ■'^t(j)) ^'^ interpreted as the trace of the product of the matrices taken in the 
same order in which consecutive integers appear in the cycle c. 
For example, tr^(xiX2) = P(j,t(xi, X2) with 



^ ^ , t(j) = (3 + (-l)-')/2, j = l,2,3,4. (2.7) 

• • • • 




Clearly, such representation is not unique, as we can realize the same polynomial po-^i by combining an appropriate 
mapping t with any element 7 G JF„ such that 7r(cr) and 7r(7) are in the same conjugacy class. The generic 



construction is to start with a partition A h n, choose (T\ G according to (|2.2p . and then define a suitable 
function t. This generic construction indicates that the monomials p^^t are essentially the monomials that appear 
in O Theorem 2], with h = I. 

Identifying x+ with x and x^ with the transpose x' of a square matrix x, define 

g^,i(xi,X2,...,x,)- J] tr j Jlx-;^'^) I . (2.8) 

cec(7) yec 

Of course, if matrix variables Xi, . . . ,Xs are symmetric, then (7^^t(xi,X2, . . . ,Xs) = p-y_t(xi,X2, . . . ,Xs). Trivially, 
for general matrix variables the polynomials coincide when 7 S ■ 
With the above notation, we state the main result of this section. 

Theorem 2.9. Let Bi, . . . ,Bs G Mmxm and let Si, . . . , Ss be N x N positive definite matrices. If Wi G 
W(Si,Si) , W2 G W(E2,-B2) , . . . , W, G W{^s,Bs) are independent and a G , then 

E(p„,t(Wi,W2,...,W,))= J2 <Z7,*(^i>---:5^K®7.t(-.7)(Si'S2,...,S.). (2.9) 

We postpone the proof to Section [51 In Section 12.21 we give a noncommutative generalization. This section 
concludes with examples, remarks and corollaries. 

2.1.4 Example of mean and variance calculation 

The purpose of this section is to illustrate ()2.9p and also to show how the coloring t{a, 7) works in a simple case. 
Example 2.10. Under the assumptions of Theorem \2.9\ formula (j2.9p gives 

E (tr(WiW2 . . . W,)) = tr(Bi)tr(B2) . . . tr(B,)tr(SiI]2 . . . S,). 

Indeed, with a = {(1, —2), (2, —3), (s — 1, — s), (s, —1)} and t{j) ~ j , we get J-s{t) = {S}. Since a @ S ~ a and 
t{a, S) = t, we get 

P<T®7,t(<T,7)(xi, . . . ,Xs) =p^,t(xi, . . . ,Xs) = tr(xi . . .Xs), 

and 

qs,ti^i, . . . ,Xs) = tr(x'i) . . .tr(x'J = tr(xi) . . .tr(xs). 

Example 2.11. Suppose Wi,W2 are as in Theorem \2.9l From Example \2.10\ we know that E(tr(WiW2)) = 
tr(_Bi)tr(_B2)tr(I]iI]2). To compute the variance o/tr(WiW2) we therefore compute E(tr'^(WiW2)) using p.9p 
for n = A with a,t as defined in ()2.7p . Nine possible 7 G ^^{t) and their contributions are listed in Table\^ 
Thus 

Var(WiW2) = (tr(SiSi)tr(B2B^) +tr(B2)tr(B2'))tr2(SiS2) 

+ (tr(BiB;)tr(B2) + tT{Bl)tT{B2B'^) + iY^{Bi)tT{Bl) + tr2(Bi)tr(S2B^) 

+ tr(B2)tr2(B2) + tr(BiBi)tr2(B2))tr((SiS2)'). 



With Bj = 



hi, 




and = /^r this gives 

Var(WiW2) = 2MiM2N^ + 2A'/i 1/2(^/1 + M2 + l)N. 



(See also Example\4.8\) 



2.1.5 Remarks 

Remark 2.12. Wishart expansions can also be modeled by planar algebras, see [19j . 



Table 1: Terms in ()2.9p for calculating E(tr^(WiW2)); here a is given by (|2.7|) . The second column is an 
illustration for Section 12.21 



7 


cr(7) 




7r(7) 


7r(CT ® 7) 


contribution 


• • • • 





• • • • 


(1)(2)(3)(4) 


(1,2)(3,4) 


tr2(Bi)tr2(S2)tr2(S]iE2) 




1 




(1)(2,4)(3) 


(1,2,4,3) 


tr2(Bi)tr(B2B0tr((I]iE2)') 









(1)(2,4)(3) 


(1,2,3,4) 


tr2(i3i)tr(i?|)tr((l]il]2)2) 




1 




(1,3)(2)(4) 


(1,3,4,2) 


tr(BiBi)tr2(B2)tr((SiE2)') 




6 




(1,3)(2,4) 


(1,3)(2,4) 


tr(BiBi)tr(B2S^)tr2(EiS2) 




5 




(1,3)(2,4) 


(1,3,2,4) 


tr(BiSi)tr(B|)tr((SiS2)2) 









(1,3)(2)(4) 


(1,4,3,2) 


tr(i?2)tr2(B2)tr((S]iI]2)2) 


;>>;:; 


5 


':XX. 


(1,3)(2,4) 


(1,4,2,3) 


tr(B2)tr(B2B^)tr((SiS2)2) 


• • • • 


4 


• • • • 


(1,3)(2,4) 


(1,4)(2,3) 


tr(B2)tr(B|)tr2(EiS2) 



Remark 2.13. Applying (jl.ip to complex Gaussian X, one can define the corresponding "generalized complex 
Wishart" matrices with complex-valued B and Hermitian E. 9, Theorem 2] can be generalized to this setting, 
and the resulting formula is 

E(p,.t(Wi,W2,...,W,))= P^-^ABi,---,Bs)p.@y.t{^i,^2,...,^s). (2.10) 

Moreover, since 7, cr G we can identify them with the pairs of permutations 7r(7) and 7r(cr) respectively; this 
identification does not affect monomials p-y which depend only on the cycle type of permutation 7r(7). When 
£ = (+, ...,+). we have ~ P-y, and t(cr, 7) ~ t; since the identification maps Brauer product to composition 
of permutations, see (|2.5p . we can write this formula directly in terms of permutations: 

E(p.,t(Wi,W2,...,W,)) = J2 P7-i,t(5i,---,S,)p„o7,t(Ei,E2,...,S.). 

Proof. Indeed, in the complex case the only pairings 7 that contribute to Wick formula are 7 G J-^ , so the 
proof of (|2.9p gives 

E(p,,t(Wi,W2,...,W,)) = J2 97.t(-Bl,...,S;)p,@^,t(Si,S2,...,S.). 



Then we note that with * denoting the conjugate transpose, we have ti{AB...C) ~ tT{C* . . . B*A*), so if 
c e C(7) then 



tr 



iiB'tu) =ti- n (V =tn n Btu) n ^^o) 



2.1.6 Special cases of interest 



In this section we list special cases of Theorem 12.91 Several similar one-matrix cases appeared in the literature, 
but without the use of the Brauer product the formulas may look differently at first. 

The following specialization of Theorem 12.91 shows that in the real-Wishart case the Brauer plays a role 
analogous to the composition of permutations for the complex- Wishart case in [51 Theorem 2] . 

Corollary 2.14. // Wi e W (Ei, /a/J , Wa G W (Ez, ^a/J , • • • , G W (S,, /a/J are independent and a e 
J-^ , then 



E(p,.i(Wi,W2,...,W,))= n^'^r'^^W.t(-.7)(Si'S2,...,E,). (2.11) 

Proof. Let M = max{Mi, . . . ,MJ. We apply Theorem [M] to M x M matrices 



hi, 




Since the cycles of Tr{-f) preserve t, a cycle of color r contributes one factor tr(_Br) = Mr to q^^t{Bi, . . . ,Bs). 
Thus q-,4Bi, . . . , S,) = n;=i M*^'^^\ 

Remark 2.15. For a single real Wishart matrix, Graczyk, Letac & Massam |17j compute moments of more 
general monomials that involve additional arbitrary matrix parameters. It is plausible that Corollarv 12. 141 can 
be deduced from their result by consecutive integration with respect to some of the auxiliary matrix parameters. 

The following one-matrix case of Theorem 12. 91 is a version of [IHl Theorem 3.5], who did not use the Brauer 
product explicitly. (In the one-matrix case t = 1 and q^^t = Py.t, so we will write Pj for both polynomials.) 

Corollary 2.16. //W G W (E,i3), B is symmetric, and a G J-^ , then 

E(p,(W))= P,iB)Pa@y{^). (2.12) 

To compute p^{In) one only needs to count the number of cycles in 7 G Tn- So when S, B are both identity 
matrices, the expression depends only on the number of cycles of the appropriate 2-regular graphs. 

Corollary 2.17 Corollary 2.4]). For A (Ai, A2, . . . ) h ?i, let pa(x) = tr(x^i )tr(x^2) . . . . // 

W G W {In,Im), and ax G is given by (|2.2p . then 

E(pa(W))= Y M#'^(*'^^)iV#'^('^^"-). (2.13) 

Proof. From the definition of 7r(7) it is clear that the number of cycles of (5 U 7 matches the number of 
cycles of 77(7). 

The cycles of 7r(o' ® 7) arise from adding the vertical edges to connect the first and last row of the three-row 
diagram. In view of our choice of a\, this gives the same cycles as adding to the graph of 7 the edges from the 
graph of (T defined by axii) — — o'>(— i), which is just the graph of ax with swapped rows. Of course nlax) is the 
inverse of TT{ax) so it has the same cycle structure. Thus px — Pax- Applying (|2.9p to Pax^ end with (|2.13p . 



2.2 Noncommutative version 

To extend Theorem 12.91 to noncommutative setting we need additional notation. 

For 7 G Tn, let cr(7) denote the number of crossings of 7 treated as the pair partition of the ordered set 
(1, —1, 2, —2, . . . ,n, —n). (The second column of Table [T] lists some examples.) 

In order to extend the definition of polynomials Pu^t to matrix arguments with entries that do not commute 
we need to specify the order in which the traces appear. To do so, we follow the standard conventions: 

• wc write the cycles of permutation T:{a) starting with the smallest clement. Thus if Ck = {jkiijki, ■ ■ ■ ijkik) 
is the fc-th cycle, then we assume that jki < j for all j G Cfc. 

• we order the cycles ci, C2, . . . by the order of their first elements. Thus 1 = jn < jai < . . . . 



Under these conventions, we define 

p.,,(W:, . . . , W,) = tr(W,(i)W,(,.,,)W,(,.^,^))tr(W,(,,^)W,(,,,) . . . W,(,,,^)) . . . 
We will use shorthand notation 

p,,t(Wi,...,W,)= n tr|nW,(,) 

cGC(a) yec 

which should be interpreted in the above sense. 

Our next result extends Theorein l2.9l to noncommutative setting. We remark that due to noncommutativity, 
we allow only special choices of ct, and that polynomials q-^^t are not needed, as matrices Bi, . . . , Bs are symmetric. 

Theorem 2.18. Let Bi, . . . ,Bs G Mmxm{^) and Si, . . . , Ss G A^jvxiv(R) be positive definite matrices. Let 
<7eT+ be given by ^ and lett : {1,...,??} ^ //Wi G (Si,Bi),W2 G >v/(E2, S2) , . . . , W, G 

W(j {Y,s,Bs) are q-orthogonal, —1 < q <1, then 

T(p,,i(Wi,...,W,)) = Yl <Z"^^^P7,t(Sl---^«K®7,*(-,7)(Sl'S2,...,S.). (2.14) 

7e^™(t) 

The proof is given in Section [31 Here we apply (|2.14p to deduce a non-asymptotic representation of the 
Marchcnko-Pastur law. 



2.2.1 A matrix representation of the compound Marchenko-Pastur law 

There is a number of results that relate Marchenko-Pastur law to the semicircle law. The well known result 
from [37] says that if X has the semicircle law then has the Marchenko-Pastur law with parameter A = 1. 
(In fact, Ref. [37] shows that for every value of A > 0, one can find a projection P such that XPX has the 
Marchenko-Pastur law.) A related representation arises by interpreting an nM x nN matrix of independent 
random variables as an M x TV block matrix; the limit as n ^ cw represents the Marchenko-Pastur law with 
parameter M/N as X*X with a matrix X of free circular entries. Yet another representation that uses free 
circular elements appears in [431 Proposition 5.2]. 

Our representation of the Marchenko-Pastur law differs slightly, as our matrix X is not self-adjoint and has 
semicircular entries. Of course, the result is not surprising; but an argument is needed as occasionally squares 
of symmetric and nonsymmetric matrices behave differently, see |34[ Section 4] . 

We first recall some properties of non-crossing partitions and recall the definition of a compound Marchenko- 
Pastur law. 

A partition V of an ordered set is non-crossing if cr(V) = 0. By NC[1, . . . ,n] we denote the set of all 
non-crossing partitions of (1, . . . , n) . 

The free cumulants {rk '■ k > 1} of a probability measure fi with all moments are defined implicitly by 
requiring that for all n G N, 

/ x-fi{dx) = E n 



veNcli,...,n] vev 



For more details, see [21] Section 2.5]; for the operator- valued case, see [41] 
The Marchenko-Pastur law with parameter A > is the law 



/,N / ^^+^ v/4A - (x - 1 - A)2 

TTxidx) = (1 - A) + do + — ^il-^)2<a,<(l + ^)2dx, 

where i5o denotes the point mass at 0. The free cumulants of the Marchenko-Pastur law are = A for all fc > 1. 

If J/ is a probability measure on the Borel sets of M with finite moments, then the compound Marchenko- 
Pastur law Tr^,x{dx) is defined as the law with free cumulants 



rfe = A / X ii{dx). 
Thus the n-tli moment of the compound Marchenko-Pastur law is 

/ x"7r,,A(dx) = V n / x*''^Iidx), (2.15) 



veAfC[i...n] vev 



compare [HJ (3.3.10)]. In particular, the Marchenko-Pastur law corresponds to = 61. 



Proposition 2.19. Let v he the distribution of the eigenvalues of a real positive- definite matrix B with respect 
to the normalized tiM '■ -MmxM — > K- (That is, v assigns mass 1/M to each of the eigenvalues of B). Let 
A = M/N. If "W G Wq {In,B) and q = then has compound Marchenko-Pastur law ■Kv,\ with respect to 

state T o tr^v- In particular, i/W S Wq {In,Im) and q = then -^W has Marchenko-Pastur law with parameter 
A. 

A pair partition 7 G J-n is non-crossing, if cr(7) ~ 0. Notice that all non-crossing 7 are in J-^ . Indeed, our 
ordering of {1, —1, . . . , rt, —n} implies that any two numbers of the same sign have an odd number of elements 
between them. Thus any pair in 7 of the same sign must generate a crossing. 

We will need to know that for a non-crossing 7 G .F„, permutation 7r(7) is the inverse of a non-crossing 
permutation. Biane [3] defines non-crossing permutations f3 G Sn by the following two requirements: 

(a) Each cycle of /3 can be written in increasing order 

(b) The partition of {1, . . . , n} into the cycles of /? is non-crossing. 

Biane characterizes non-crossing permutations by the geodesic condition 

#C(/3)+#C(po/?-i)=n + l. (2.16) 

(See [33 Remark 2.9], [HI Section 6], and [H Exercise 6.19.hh]. Note that this is g in (jO]) below.) 

It is well known that the non-crossing pair partitions of 2n elements are in a one-to-one correspondence 
with the non-crossing partitions of n elements, see [42l Exercise 6.19]; see also [211 page 74]. It turns out that 
the restriction of n : J-n — > Sn to non-crossing pair partitions defines an explicit version of this bijcction. 

Lemma 2.20. Let tt be the mapping from pair partitions of {±1, . . . , ±n} to partitions of {1, ... ,71} which 
assigns to j ^ J-n the partition into the cycles of permutation tt{j). Then n is a bijcction from the non-crossing 
pair partitions with respect to ordering (1, —1, . . . , ?i, — ?i) to the non-crossing partitions 0/ (1, ... , n) in natural 
order. Moreover, for non-crossing 7 G J^n, 0,11 cycles of Tri^) can be represented as decreasing sequences. 

(We note that our "standard" representation of a cycle begins with the smallest element, which therefore 
needs to be moved to the end of the cycle to get a decreasing sequence representation.) 
The proof below is due to the referee. 

Proof. Non-crossing pair partitions are in Tn , and we already noted in Remark 12.31 that tt : J-^ Sn is 
one-to-one. 

Let 7 G Tn and let a G 5 be such that 7r(7) = a. Notice that permutation 7 o (1, —1, 2, —2, . . . , n, —n) works 
by sending i ^ a~^{i) and —i —(a o (1, 2, . . . , n)), for every 1 < i < n. Hence 

#C (7 o (1, -1, 2, -2, ... , n, -n)) = #C {a'^) + #C (a o (1, 2, . . . , n)) . (2.17) 

By using this fact and Biane's geodesic condition (|2.16p . we see that 7 is non-crossing with respect to 
(l,-2,2,-2,...,n,-n) iff 

#C(7) + #C (7 o (1, -1,2,-2,..., n, -n)) = 2n + l. 

Since #C(7) = n, from (|2.17p we see that this is equivalent to 

#C (a-i) -f #C (ao(l,2,...,n)) =n + l, 

which, again by (|2.16p . is equivalent to the statement that a^^ is non-crossing. 

Proof of Provosition \2.1S\ When g = 0, (|2.14p gives 

iV-"r o tr7v(W") ^ ^ ^#c(<T®7)-n-i "Q tr(S#"), 
7 cec(7) 

and the sum is taken over the non-crossing pair partitions 7. Non-crossing 7 are in J-n , so we can use p.Sp . 
With p = 7r(<T) = (1, . . . , n), we get 

Af-"rotrAr(W") ^ ^#C(po^(7))-n-i -Q tr(B#=). 

7ej^„ non-crossing cgc(7) 



By Lemma l2.20| we can take the sum over the permutations a obtained as images of non-crossing 7. 



a c£C{a) a cGC(q) 

If the cycles of a G Sn form a non-crossing partition and are in decreasing order, then the cycles of (3 a^^ 
are in increasing order and form a non-crossing partition, too. Since a and have the same number of cycles, 
(|2.16p can be written as 

#C(poa) + #C(a) =n + l. (2.18) 

Thus, 

Qe5„,Q-i non-crossing cec{a)-^ veNCli,...,n] vev^ 

which matches the moments p.lSp of the compound Marchcnko-Pastur law. 



3 Proofs of Theorems [HI and [27181 

Following [33| . see Lemma 1], we write polynomials p^^t and q-y^t as follows. 
Lemma 3.1. Let 7 e JF„. With a = n{-f), (ei, . . . e„) = £(7), 

n 

p^,t(xi, . . . ,X,) ^ Y[[xt{t)]l{i),I{a{i)) , (3.1) 

I:{l,...,n}^{l,...,N} i=l 

n 

/:{l,...,ri}-f{l,...,Af} 1=1 

Proof of Theorem \2.9l Denote p = 7r(cr) G 5„. Expanding matrix multiplications in (j3.ip . we write the 
polynomial on the left hand side of (|2.9p as 



p,,t(Wi,W2,...,W,) = ^ [][4(,)X;(,)St(,)X4(,)At(,)],(,),^(p(,;)) 

/:{l,...,n}-»{l,...,JV} i=l 

(n \ ^ ^ 

1=1 J i=l i=l 

where the sum is now taken over all functions I ,U,V : {1, . . . , rt} ^ {1, . . . , and over all functions J : 
{1, . . . , n} ^ {1, . . . , M}. The Wick formula expresses the expectation of the product of jointly normal centered 
random variables as the sum over all pairings of product of the covariances for each pair of random variables, 
compare Definition 11.21 with q = I. To describe the pairings, we associate positive integers j = 1,2, ... ,n with 
'^t{j) ^'^'^ negative integers —j with Xj(j). Then the pairings are determined by 7 G !Fn{t)] denoting by 
a = 7r(7), Table [D lists four types of matches that might contribute to the moments. 
The expectation becomes 

^\X\\^'t(i)\u(i).K{i)\y^t{i)]j(z)y(i)\ = Y Xa{'1,U,V)xb{j,K,J). 

where Xa(7i f^i ^)X-b(7i^j J) = 1 only when functions U, V and K, J, fulfill the matches specified in the last 
column of Table [2] at each edge of 7; otherwise, xa{'^, U, V)xb{j, K, J) = 0. 
This allows us to write 

E(p,,t(Wi,W2,...,w,))= n^(7)ni3(7), 



Table 2: Types of matches along 7 S J^n{t); here a = 7r(7) G 5„ 



Match type 


graph 


matching entries 


identification of indices 


A 


• • 


i'^'tU)]uU)KU) = [^t(a{j))]j{aU))V{a{j)) 


U{j) = V{a{j)), K{j) = J{a{j)) 


B 


• • 




U{j) = U{a{j)), K{j) = K{aij)) 


C 


• • 


i^tU)]jU)VU) [^t{aU))]j{a{j))V{a{j)) 


J{j) = J(aO-)), V{j) = V{a{j)) 




• • 

• • 


\^tU)].JU)V{]) [^'t{a{j})h(aij))K(a(j)) 


J{j) = A^(a(j)), = f/(aO-)) 



where 



Ha (7) := 5^ XA(7,C/,^^)n[^^w]^(0.c/w[^ 



t(»)JvW,/(p(i)), 



It remains to prove the following two results. 
Lemma 3.2. If j ^ ^n{t) then 

11^(7) = P<T®7,t(,T,7)(Sl, . . . , Ss). 

Lemma 3.3. 7/7 e t/ien 113(7) — Qt,t{Bi, . . . ,Bs). 



Proof of Lemma [ 

Given K, J : {1, . . . , ,n} {1, . . . , M}, define K : {±1, . . . , ±n} ->{!,..., M} by 



ifj>0, 
J{-j) ifj<0. 



It is clear that {K, J) t— » is a bijection between the set of all pairs of functions K, J : {1, . . . , ,n} ^ {1, . 
and all functions K : {±1, . . . , ±n} ^ {1, . . . , M}. 
From Table [H it follows that 

XBh,J,K)= Y[ 

{u,v}£f 

For a fixed 7 £ denote a = tt{j), (ej) = £(7). Then we can factor 115(7) over the cycles of a, 

n 

^B{l)^^XB{l,K,J)W[Bt^,)\Ki^),J[^)= H (7, ^ I c ' -^D H ^^O^^U) ' 

K,J i=l c<£C{j} K\^,J\^ jec 

On a single cycle c = (ji, . . . ,jr) of a, t := is constant. So to end the proof, it suffices to show that 



^ XBil,K,J)Y[[Bt]KU)ju)=i'^ 

K,J:c^{l,...,M} j&c \j€c 

We note that cycle c corresponds to the sequence 

C = (jl, ±j2), (TJ2, ±J3), ■ • ■ , (Tjr-1, ±jr), (Tjr, -jl) 



of the ordered edges of 7. The corresponding edges of 6 are 

S\c = ((-il,jl),(±j2,Ti2),-.-,(±jr,Tjr))- 

These two sequences can be written directly in terms of (e-,) as 

C = (ji, -ei2J2), (ej2J2, -EjaJs), ■ • ■ , {£j,_Jr-i, (e>>, -ji)- 

For fixed K, J, the expression under the sum on the left hand side of (j3.7p can be written as 



1=1 

r 

(Here we use cyclic convention jr+i ;= ji-) Thus the left hand side of p.7p is 



K:{±J1,...,±>}-{1,...,A/}' 



..,M} »=1 

proving (|3.7p . 



Proof of Lemma [ 

We begin by rewriting 11^1(7) in terms of matrices Si, . . . , Sg. Recall that p — n{a) £ 5„. We claim that 
with I: {±1, . . . , ±71} {1,...,N} defined by 

and with t : {±1, . . . , ±n} — > {1, . . . , s} extended by t{u) :— t{\u\), we have 

nA(7)= E n Mnuui.r (S-H) 

/:{!,. ..,n}-»{l,...,Ar} {u,t>}G7 

To prove p. lip , we re-paramcterize each pair of functions U,V by a single function L/ : {±1, . . . , in} 
{l,...,iV} defined by 

By inspecting the rows of Table [H we see that 

{M,u}e7 

For fixed /, 

n 

U,V i=l 



factors over the cycles of a into the product of 

H XA(7,C^,^)n[^t(j)]^(i)'C^O')[^tO)]^0)'^(pO))- 

lJ,V:c^{\,...,N} jGc 

The corresponding set C of ordered edges of 7 is 

C = (Ul, Wl), (m2, W2), . . . , (Wr, Wr), 

and, as in the proof of Lemma 13.31 we have Ufe+i = — Wfe, /s = 1, . . . ,r with it^+i := ui > 0. Of course, 
c = (ji, . . . = (I^^iIj • ■ • j l^^rl), and f is constant on c. Thus, using the symmetry of 

r 

jGc fe=l 



fc=i 

Since = 1 imphcs U{uk) = U{vk), we see that 

[/(wfe+i) = U{-~Vk). 

Thus, 

E xA(7,c/,^)n[^t (j)]^0).c/0")[^t(j)]v'(j),/(p(i)) 

C/,y:c-f{l,...,Ar} j'Gc 

r 

X] \\'^U{vk+i) = U{-Vk)^Ml(^-vk),U{-Vk)'^Mu{vk+i)J{vk+i) 

U:{±Ju...,±jr}^{l.....N} 

H [^*]7(-i>,),c7(-t,fc)[^*]c7(-t,fc)JK+i) = n P*]/(«),7(t.)- (3.14) 

;7:{-ui,...,-ur}^{i,...,Af} {u,t)}ec 

We now fix 7 £ !Fn{t) and consider p. lip , which we want to factor over the cycles of /3 = ■7r((T ® 7). To do 
so, let r be the bijcction from pairs {u, w} G 7 onto the corresponding pairs {r^s] £ cr @ 7 which appeared in 
Definition [Ml see Table O 

Since r is a bijection and the cycles of (3 are in one-to-one correspondence with the sequences of pairs in 
cr ® 7, to prove (|3.4p it suffices to show that if c is a cycle of /3 and C is the corresponding sequence of pairs in 
cr @ 7, then 

n P*h]7(«),7(«) = W\^t{<yn)U)\iU)J°m (3-15) 
{M,t.}er-i(c) jec 

We prove (|3.15p by showing that each of the factors on the left hand side of (|3.15p can be rewritten into 
the corresponding factor on the right hand side of p.lSp . To do so, let {w, u} = r~^({r, s}), where {r, s} £ C 
are such that |r|, \s\ arc the consecutive numbers in cycle c, i.e., \s\ = /3(|r|). 

We note that r preserves the number of positive elements in each pair {u, w}; thus if u, w > then r, s > 
and similarly if < then rs < 0; Table [3] lists all possible cases with the intermediate three-row diagram. 

Wc now consider all possible cases. 

A-1 Suppose r < 0, s > 0. Since T,t{u) = ^t{v) is symmetric, without loss of generality we may assume 
u > 0,v < 0. Then u = s and r = —p{—v), so p.lOp gives 

Pt(")]7(«),7(«) = Pt(-")]7(«),7(«) = Piop-i(k|)]/(|r|)J(s) = [^t{a,j){\r\)]l{\r\)Jc.l3{\r\]- 

A-2 Suppose r > 0, s < 0. Since St(u) = ^t(u) is symmetric, without loss of generality we may assume 
u > 0,v < 0. Then u ~ r and s = —p{~v), so p.lOp gives 



P*(")]7(u)J(t,) = Pt(|r|)]/(r)J(|s|) = Pt(^,7)(|r|)]/(|r|)Jo/3(|r|)- 



Table 3: Bijection r from pairs {u, t;} G 7 to pairs {r, s} S cr ® 7. Here p = 7r(cr) € 5, 



Case 


{u,v} e 7 


{r, ,s} = t{{u,v}) 


A 


« ... . p(-f) 

* ... • • 


{r, s} = {li, -p(-u)} 


B 




{r, s} = {w, u} 




— ^ • . . . -V , 
• • 




C 


u . — — ■ ' ■ , V 

• -p{-u) ... . -p{-v) 
• • 


{r,s} = {-p(-m), -p(-w)} 



B Suppose r > 0, s > 0. Then {w, w} = {r, s} and (|3.10|1 gives 

Pt(f)]7(u),7(t,) === Pt(r)]/(r),/(s) = [St(|r|)]/(|r|),7o/3(|r|), 

C Suppose r < 0, s < 0. Then swapping the roles of u, v if necessary, we may assume that r = —p[—u) and 
s — -~p{—v), so p.lOp gives 

\-^tiu)]T{u)J{v) = Ptop-i(|r|)]/(|r|)J(|.s|) = [^t{a,f){\r\)]l(,\r\)JoP{\r\}- 

Proof of Theorem \2.18[ This proof is similar to the proof of Theorem 12.91 with p.2p used instead of Wick 
formula. Throughout the proof, we denote p ~ 7r(cr), P = 7r(cr @ 7). For functions / : {1, . . . , n} {1, . . . , A^} 
and J : {1, . . . , n} {1, . . . , M}, we denote by /, J functions on {±1, . . . , ±n}, which are defined by ()3.10p and 

■— '^(IjI)i respectively. We also extend t to {±1, . . . , ±n} by t{s) :— t{\s\); we will use the same symbol t 
for this extension. From (|2.2p and (|1.3p we see that 

T (p^,t(Wi, . . . , W,)) = ([^*(l)]j(l),7(l)[^*(l)b(-l)J(-l) • • • [^t(«)]j(«)J(n)[^t(«)],7(-n),7(-n)) ' 

By pT^ . this implies 

Since ^[^t(ti)] [^((■u)] 7(„) 7(^,)^ =0 when t(it) ^t{v), we see that the sum can be restricted to the 

color-preserving partitions 7 G J-n{t)- Changing the order of summation, and taking into account (|1.2p . we get 

r(p.,(Wi,...,w,))= 9"^"^E n E 9"-(-'ni(7)n2(7), 

where 

ni(7)= E n [^*H]j(n).7(.)' (3-16) 

J:{l,...,n}->{l,...,j\/} {u,v}e-f 

"2(7)- E n P*(«)]7(«).7h- (3-17) 

I:{l....,n}^{l,...,N} {u,v}e-f 

To end the proof, it remains to show that 

ni(7)=P7,t(-Bi.----5.0 (3-18) 

and 

n2(7) = P<T®7,t(<T,7)(Sl, . . . , Ss). (3.19) 

The proof of (|3.18p is similar to the proof of Lemma [3.31 and is omitted. Since the right hand side of p.l7p is 
the same as the right hand side of p. lip , formula (|3.19p follows from the proof of Lemma [321 



4 Proof of Theorem 11.71 



The identification of the limit in Theorem 11.71 when q = or q = I is a consequence of the following general 
property of g-Wishart matrices. 

Theorem 4.1. Let 

wr ^ , wf ) , . . . , w(^) , wr \ ) , . . . , w(^) e ( 1/^ , /m) 

be q-orthogonal with M ~ M{N) such that M/N ^ X. Fix a polynomial Q{xi^...,Xs) in noncommutative 
variables and let 

^ tr(Q(W(^\ . . . , WW)) - r(tr(Q(w(^', . . . , W^))) 

and 

= tr(Q(wl^\ . . . , wW)) - r(tr(Q(w(^\ . . . , Wf)))). 

Then for m > the limits limAr_,oo t{X'^), limjv-»oo t H^n + Yn)"^) exist, and 

hm (r {{Xn - Yn)HXn + Y^D - 2t{XI)t {{Xn + Y^D) = 0- (4.1) 

Remark 4.2. Since the conditional expectation of X^ given Xjv + i^jv is {Xm + YAr)/2, expression (|4.ip can be 
interpreted as the statement that the conditional variance of Xf^ given X^ + Y^ is asymptotically constant. In 
the self-adjoint case, the latter property identifies the distribution when g = 0, 1, see Theorem 3.1]. We note 
that Theorem 14.11 holds true for all — 1 < g < 1 , but Examples 14.61 14.81 and 14.71 show the limiting law may fail 
to be g-Gaussian for < jql < 1. This suggests that there is no natural version of [5]. 

4.1 Auxiliary combinatorial result 

The proof of Theorem 14.11 relics on the following partition version of a theorem of Jaques [22] . 

Proposition 4.3. // 7 G Tn and a G J-^ are such that the graph 7 U (t U (5 is connected, then there exists an 
integer h >0 such that 

#C(7) + #C(a) + #C(a @ - n = 2 - h. (4.2) 

For the proof, we will need to consider J-n in all three roles: as pair partitions, as graphs, and as a subset of 
S2n- To help distinguish between these interpretations, denote by z{a) the number of cycles of a permutation a. 
Recall that in this notation, if 7 £ JF„ then #C(7) = z{tt{j)). We will use the above mentioned result of Jaques, 
which we now state in the version most convenient for our purposes. 

Theorem A. If a, P £ iS„ act transitively on {1, ... ,n}, then there exists an integer g > such that 

z{a) + z{P) + z{p-^a) -n = 2-2g. (4.3) 

(For an accessible algebraic proof see j29| Theorem 1]. For a graph-theoretic proof, see [Ul Proposition 
1.5.3].) 

Proof of Proposition \4-3\ The proof consists of two claims, and an application of ()4.3p . As noted in the 
proof of [ini Lemma 3.5], 

Claim 1. For 'j £ Tn, 

z(7o5) = 2#C(7). (4.4) 

(This is because every cycle of (5 U 7 splits into two cycles of 7 o 5. Recall that S is defined by ()2.ip .) 
We now note the following. 

Claim 2. If j G J-n, c G J'n '^'^'^ graph 7 U cr U (5 is connected, then the group generated by permutations 
J o S,a o S has at most two orbits on {±1, . . . , ±n}. 

To prove this claim, we note that for a G J-^ the cycles a o 5 arc exactly the cycles ci, . . . , Cfc of 7r(cr) G Sn 
and — ci, . . . , — Cfe. 

Since the graph 7 U ct U J is connected, 7 must connect the sets 



Ci U -Ci, C2 U -C2, . . . , Cfc U -Cfc , 



so "f o 5 must also connect these sets. Thus if S is an orbit of {a o S,j o S), then {|j| : j G 5} = {1, . . . , ??}. 

As the mapping S 3 j ^ \j\ G {1, . . . , n} is two- valued and onto, there cannot be more than two disjoint 
inverse images of {1, . . . , n}. 

By ClaimO applying (|4.3p toa — joS,(3 = aoS£ S2n or to their restrictions on the two orbits separately, 
we have 

r^ , // s:\-i x\ o (2-2.9 if (cr o J, 7 o 5) acts transitively 

I 4 — 2gi — 2g2 otherwise 

for some integers 5,(71,(72 > 0. Thus in both cases, there is an integer h>0 such that 

z(7 06) + z{a o (5) + z((cr o o 7 o (5) - 2n = 4 - 2/i. (4.5) 

We now note that (ct o 5)^^ ~ & o 5, where (t(j) = Moreover, (TO(5 = (5o(tso permutation = (<t o 5)^^ o 

7o^ = i7o5o70(5is similar to it o 7. 

In the proof of CoroUarv 1 2 . 1 71 we already noted that #C((T @ 7) is the number of components of the graph 
(J U 7. We now observe that the cycles of (i o 7 split each such component in two. Therefore 

z{{<7 o (5)^-'- o 7 o (5) — z[g 07) = 2#C(ct ® 7). 

This, together with Claim [T] and (|4.5p gives ()4.2p . 



4.2 Proofs of Theorems HI] and O 

Our first goal is to show that the limits of joint moments exist; the main step consists of analyzing moments of 
the monomials. 

We will need additional notation. By JF„(t, a) we denote the set of all 7 G Tn^t) which connect every cycle 
of a to some other cycle. That is, all components of the graph 5 U (T U 7 are strictly larger than the components 
of cr U 5. 

Lemma 4.4. Let Wj^, W2, . . . , G Wg (-^/at, Im^ be q-orthogonal. Suppose that the cycles Cu ^ u < r) of 
cr G J-^ consist of consecutive integers, as in (j2.2p . Fix t : {1, . . . , ?i} —>■ {1, . . . , s}. Then 

r I f[ Li n W,(,)) - r(tr( [] W,(,))) ) ) = ^ ^„(,)^^#C(.)^#C(.®.)-„, (4 g) 

Proof. Under the standard convention that empty product is identity, the left hand side of (|4.6p is 

E nMiiw^o))) n^t^diw^o))). (4.7) 

(Here, finite subsets F C {1, . . . ,r} and their complements in {1, . . . ,r} are treated as ordered sets with the 
natural order on integers, so that the noncommutative products of traces that occur in the first factor are 
well defined. That is, if _F = {ui, U2, . . . , it,„} with ui < M2 < • • ■ < Um then Jl^^gj^ Xu '■— Xu-^Xu2 ■ ■ ■ Xu,„ is well 
defined.) 

We now re-rewrite the right hand side of (j4.6p . Let ci, . . . , be the cycles of a in increasing order, and let 
fci, . . . ,kr denote their lengths. (Recall our convention that ci, . . . , are in fact the cycles of Tr{cr) G Sn-) 

For 1 < u < r, let Qu denote the set of 7 G Tn{t) which "isolate" cycle c„, i.e., do not connect cycle Cu to 
any other cycle. Then Tn(t, a) = J-n{t) \ {J^^i Gu, so for any weighted counting measure H, 

r 

H{Tn{t, a)) = H{Tnit)) -H{\Jg^)= E i-^)*^H{ f| ^0, 

"=1 FC{l,...,r} u£F 

where H{!Fn{t)) corresponds to = in the expansion. Therefore, the right hand side of (|4.6p is 

Fc{l,...,r} 7en„eFe„ 



Note that 7 G HueF^" '^^^ uniquely split into the disjoint union 70 U[J^^p ^u, where 7u G J-k^ is a pair 
partition of cycle c^, and 70 G with rio = n — X^ugf ^ti- W'-' ^'^^ note that a @ 7 = (ctq ® 70) U IJueF c« ® 7m- 
Since the cycles of a consist of consecutive segments, no additional crossings of 7 can arise from the edges 
connecting Cu besides the crossings already accounted for as crossings of 7u. Thus cr(7) = cr(7o) + J2ugf '^i'(7m)i 
and the right hand side of (|4.6p can be written as 

Applying ()2.14p to each of the factors, we arrive at (|4.7p . thus deriving the left hand side of (|4.6p . 

Next, we show that the non-negligible terms in expansion (|4.6p correspond to 7 that connect pairs of cycles 

of cr. 

Lemma 4.5. Fix 7 G J^n{t, o"). Then the limit 

Mnh) := lim Af #c(7)7v#c(<T@7)-n (4 8) 

JV-»oo,A//Ar-tA 

exists. Moreover, if the number of cycles of a connected by 7 is more than 2, then M„(7) = 0. Ifj connects only 
pairs c-u,c-u of cycles of a, and "fu,v denotes the corresponding part ofj, then 

u,v } ■ 

{u,v} 

(The product is taken over all unordered pairs u,v corresponding to the connected cycles of a). 

Proof. Since 7 G .F„(t, cr), once we decompose cr U 7 U (5 into connected components, the decomposition 
splits the cycles of cr into groups of mi, . . . , > 2 cycles, perhaps non-consecutive. This partitioning splits 7, cr 
and cr ® 7 into disjoint triplets cr„, cr„ ® 7„ G Tn^ , m = 1, 2, . . . , r, with 71 = 711 + • • • + n^, which contribute 
multiplicatively, 

r 

«=i 

(Recall that 77i„ = #C(cr„).) Since the graphs of cr„ U 7„ U are connected, by formula (|4.2p . there are integers 
hi, . . . hr > such that 



M#C(7)jV#C(^@7)-n ^ (A//Ar)#C(7) -Q TV 



iv 

u=l 



Therefore, the limit always exists. The limit is zero if at least one of the 7Ti„ > 2 or at least one of the hu > 0. If 
all groups arc pairs, i.e., mi = 7712 — ■ ■ ■ — = 2, then the limit factors into the contributions from the pairs 
of cycles of cr; of course, the only 7 that contribute to this limit are those with hi = h2 = ■ ■ ■ = hr = 0- 

Proof of Theorem \4-.l\ Since {x — y)'^ = x{x — y) + y{y — x), and joint moments of [Xn, Yn) are symmetric 
in X,Y, it is enough to show that 



N 



lim (r {Xn{Xn - Yn){Xn + YnD - t{XI)t{{Xm + Iat)")) = 0. (4.9) 



After re-indexing the g-Wishart matrices so that W_j W-,-, we can represent Xjv = "^"^a^," ^'^'i 
Y]s[ ~ ^ ctuYN,u as finite linear combinations (with the same coefficients, and the same choice of functions 
tu '■ Cu {1, . . . , s}) of the centered monomials 



Xn,u = tr( n W,„(,))-T(tr( n W,„(,))), (4.10) 
Yn,u - tr(n W_,„(,))-T(tr(n W„,„(,))). (4.11) 



(Here c„ are finite sets of integers of prescribed cardinality, which will be assembled into the cycles of a in the 
next step.) 



Expanding f|4.9p into the monomials, it is enough to show that 

m m 

Yirn^ {t{Xn,v,{Xn,u - Yn,u) Y[{Xn,v, + - r(Xjv,„oXjv,n)r( Y[{Xn,v, + Yn,v,))) = (4.12) 

for all admissible indexes vq, . . . , v„i, u. Once wq, u, wi, . . . , Wm are fixed, we modify t and the sets c„ so that 
CvotCutCv-^^ . . . ,Cv,^ are consecutive blocks of integers that form the cycles of cr € J-^- This allows us to apply 
Lemmas 14.41 and I4.5[ which imply that the limits of moments exist. Expanding the products 

n m 

n (^A^.-. + yn..., ) - e n ^N,., (4.13) 

j = l i)i,...,?;„=0,l j=l 

we see that it suffices to show that 



^lim^ E ^ ^^^-0 {Xn,u - ym,u) n ^N.., yk'; 

i7i,...,t;„=0,1 y i=l 

= }'Zo ^ <XN,.,X^.u)r{\{x]^^^,Yl,~:;). (4.14) 

77i,...,r/„i— 0,1 j — 1 

To shorten the notation, denote by cq, c, Ci, . . . , Cm the cycles c„p, c„, c^^ , . . . , Ci,^ of a. In view of (|4.10M.l"T|) 
Lemmas [4.41 and ITSl implv that for each choice of 77 = (771, . . . , rjm) S {0, 1}" there is : {1, . . . , } ^ {±1; • ■ • i ±s} 
such that 

m 

r(X^,„„Xjv,nn^^%^J^r.?) = E '7"^^^M™+2(7) + 0(l/7V). (4.15) 

For j = 0, . . . , m denote by Ti-jitri) the set of 7 G ^^(i,,, cr) which connect cycle c to cycle Cj and do not connect 
the pair c,Cj to any other cycles of a. Lemma implies that ()4.15p can be expanded as 

m m 

r(X^,„„X^,„n^;^%,^^;'0 E 9"^"^^W2(7)+E E 'Z"^"^M™+2(7) + 0(l/iV). (4.16) 

J=l 76'Ho(t^) i=l 7e'Hi(*>,) 

All 7 G 'Ho{Uj) can be decomposed into 70 U 7', where 70 pairs the cycles co,c only, and 7' is the pairing of 
the elements underlying cycles ci, . . . , c„. Furthermore, since cq, c, ci, . . . , are consecutive blocks of integers, 
cr(7) = cr(7o) + cr(7'). Therefore q"('^)M„+2 (7) = q"^'^°'> M2i-fo)q'''^'^"> M„{j') and Lemmas|431and[13]together 
imply that the first sum simplifies, 



-VJ^ 



+ E E q"^''^Mra+2h)+0{l/N). (4.17) 

For i G {1, . . . , n}, let 

l^?7(i) otherwise. 

Since yjv.u differs from X^v^u only in the sign of t\c, and since TCo{trj) = 0, the same reasoning as above gives 

rn m 

r{XN,voYN,uY[^N,vYN7;)=T. E q"^'^M^+2h) + 0{l/N). (4.18) 

In view of (|4.17p and (|4.18p . to prove (|4.14p . it is enough to show that for every j > 1, 

E E q"^^^M,„,+2{l)=J2 E 9"'<"^M„+2(7)- (4.19) 



If 7 e Hj{t^) then (jil^ gives 



M™+2(7) 



if = 0, 



where 7j is the pairing on c U Cj and 7^ is the remaining set of pairs. On the other hand, if 7 G Ti.j{tjj) then 

M„+2(7) 



iff7i = l, 

/•/2(7,)Mm(7;) if^, =0, 



as the joint moments of X^r „, Xjv,cj are the same as the moments of Iat^uj ^at.cj and t = t outside of c U Cj. The 
later also implies that with t]' = rj except for one entry 77' (j) = 1 — rjij), we have 

Thus summing over 77' on the left hand side of (j4.19p we get the right hand side, ending the proof. 

Proof of Theorem \1.T\ If A = then all moments of Zm,n converge to zero, which we will consider as a 
normal or a semicircle law of variance 0. Through the remainder of the proof we assume that A > 0. Since 

(Q(Wi, . . . , W,) - r(Q(Wi, . . . , W,)))™ 

can be written as a linear combination of centered monomials of the form that appears on the left hand side of 
(|4.6p . combining Lemmas 14.41 and 14.51 we see that all moments converge. Furthermore, from (|2.14p we see that 
the n-th moment grows at the rate bounded by a multiple of 4i^J-n = 1 • 3 • . . . • (2n — 1), so in the self-adjoint 
case the limiting law is uniquely determined by moments. 

When q = \ and Q is real symmetric, the two copies X^r, of Zm,n are independent and converge jointly 
to a pair of independent identically distributed self-adjoint (i.e., real) random variables X,Y. From Theorem 
WJ[ we see that E {{X - Y)'^piX + Y)) = 2E(X'^)E {p(X + Y)) for all polynomials p. By the previous estimate 
for the moments, the characteristic function of X -I- y is analytic at 0, so the moment problem is unique and we 
can replace polynomial p by any bounded measurable function. Thus Var(X|X + Y) = const, which is known to 
imply that X is normal, see [26l Corollary 4.1]. 

When q = and Q is real symmetric, the two copies X^^Y^ of Zm,n are free and converge jointly to a 
pair of free identically distributed self-adjoint noncommutative random variables X,Y . From Theorem 14. 11 we 
see that r ((X — YYp{X + Y)) = 2t{X-^)t {p{X + Y)) for all polynomials p. For q ~ 0, the moments of X are 
bounded by the number of non-crossing elements in J-'„', this is the 2n-th Catalan number, so the support of 
X + Y is bounded and we deduce that Var(X|X + Y) = const, which is known to imply that X has semicircle 
law, see [5J Theorem 3.2] for a more general result. 



4.3 Concluding examples 

The following computer-assisted calculations give low order moments of some polynomials in g-Wishart matrices. 
Together, they illustrate that the limit law depends on the polynomial in a nontrivial way. 

Example 4.6. In the simplest case Wjv G Wg (-^/Ar,/j\f) and M ~ XN , we have T(tr(WAr)) = AA^. The first 
six moments for the limit of Xjy ~ tr(Wjv) — NX are consistent with the q^ -Gaussian law: with := t(X^) = 
(1 4- q)X, we have t{X^) 0, 

T{Xt,)/s' ^ (2 + q'), t{X%) ^ 0; t{X%)/s' ^ 5 + 69^ + 3g« + q^\ 

Example 4.7. To see whether one could identify the limit law for a second-order polynomial in W G 
Wg (-j^/jv, Im) , we computed the first moments of the limit as N ~> 00 and M/N —t X of 

Xn = tr(W2) - atr(W) - T(tr(W2) - atr(W)). 

The general expression for the asymptotic variance is cumbersome, but it simplifies to t{X'^) X^{q^ -\- q'^ -\- 
q^ -\- 1) when a = \ -\- q^ -\- 2A. For this value of a, the normalized fourth moments r {X%) /(''"(X^))^ converge 
to 2 + q^^ , which is the fourth moment of the q^^ -Gaussian law of unit variance. 



Example 4.8. Suppose Wi,W2 G [j^In,Im) are q-orthogonal. Then r(tr(WiW2)) = AP/N, see Exam- 
pleKM and from Table{^we see that Yai{tT(WiW 2)) = {1 + q^)q'^AP /N^ + 2{1 + q)M^ /N^ + {1 + q)q'^ /N^ . 
Denote Xn = tr(WiW2) - NX^, and suppose M = M{N) = XN + 0{1/N). Then liniAr^oo T{Xjf) = 
(g^ + q^ + 2X + 2qX) , and computer-assisted calculations give limjv_>oo T{Xff) — and 

\im^T{X%) = X-" [q^l + q^f {2 + q'^)+Aq'' (1 + q) {l + q^) {2 + q^) X + A {I + qf {2 + q^) A^) . 

After normalizing by the variance, we see that in this case the fourth moments converge to a rational function 
ofq- 
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